CRISPR Screening Uncovers a Long-Range Enhancer for ONECUT1 in Pancreatic Differentiation and Links a Diabetes Risk Variant.

Functional enhancer annotation is a valuable first step for understanding tissue-specific transcriptional regulation and prioritizing disease-associated non-coding variants for investigation. However, unbiased enhancer discovery in physiologically relevant contexts remains a major challenge. To discover regulatory elements pertinent to diabetes, we conducted a CRISPR interference (CRISPRi) screen in the human pluripotent stem cell (hPSC) pancreatic differentiation system. Among the enhancers uncovered, we focused on a long-range enhancer ∼664 kb from the ONECUT1 promoter, as coding mutations in ONECUT1 cause pancreatic hypoplasia and neonatal diabetes. Homozygous enhancer deletion in hPSCs was associated with a near-complete loss of ONECUT1 gene expression and compromised pancreatic differentiation. This enhancer contains a confidently fine-mapped type 2 diabetes (T2D) associated variant (rs528350911) which disrupts a GATA motif. Introduction of the risk variant into hPSCs revealed substantially reduced binding of key pancreatic transcription factors (GATA4, GATA6 and FOXA2) on the edited allele, accompanied by a slight reduction of ONECUT1 transcription, supporting a causal role for this risk variant in metabolic disease. This work expands our knowledge about transcriptional regulation in pancreatic development through the characterization of a long-range enhancer and highlights the utility of enhancer discovery in disease-relevant settings for understanding monogenic and complex disease.

Decoding Heterogenous Single-cell Perturbation Responses.

Understanding diverse responses of individual cells to the same perturbation is central to many biological and biomedical problems. Current methods, however, do not precisely quantify the strength of perturbation responses and, more importantly, reveal new biological insights from heterogeneity in responses. Here we introduce the perturbation-response score (PS), based on constrained quadratic optimization, to quantify diverse perturbation responses at a single-cell level. Applied to single-cell transcriptomes of large-scale genetic perturbation datasets (e.g., Perturb-seq), PS outperforms existing methods for quantifying partial gene perturbation responses. In addition, PS presents two major advances. First, PS enables large-scale, single-cell-resolution dosage analysis of perturbation, without the need to titrate perturbation strength. By analyzing the dose-response patterns of over 2,000 essential genes in Perturb-seq, we identify two distinct patterns, depending on whether a moderate reduction in their expression induces strong downstream expression alterations. Second, PS identifies intrinsic and extrinsic biological determinants of perturbation responses. We demonstrate the application of PS in contexts such as T cell stimulation, latent HIV-1 expression, and pancreatic cell differentiation. Notably, PS unveiled a previously unrecognized, cell-type-specific role of coiled-coil domain containing 6 (CCDC6) in guiding liver and pancreatic lineage decisions, where CCDC6 knockouts drive the endoderm cell differentiation towards liver lineage, rather than pancreatic lineage. The PS approach provides an innovative method for dose-to-function analysis and will enable new biological discoveries from single-cell perturbation datasets.

ChromaFold predicts the 3D contact map from single-cell chromatin accessibility.

The identification of cell-type-specific 3D chromatin interactions between regulatory elements can help to decipher gene regulation and to interpret the function of disease-associated non-coding variants. However, current chromosome conformation capture (3C) technologies are unable to resolve interactions at this resolution when only small numbers of cells are available as input. We therefore present ChromaFold, a deep learning model that predicts 3D contact maps and regulatory interactions from single-cell ATAC sequencing (scATAC-seq) data alone. ChromaFold uses pseudobulk chromatin accessibility, co-accessibility profiles across metacells, and predicted CTCF motif tracks as input features and employs a lightweight architecture to enable training on standard GPUs. Once trained on paired scATAC-seq and Hi-C data in human cell lines and tissues, ChromaFold can accurately predict both the 3D contact map and peak-level interactions across diverse human and mouse test cell types. In benchmarking against a recent deep learning method that uses bulk ATAC-seq, DNA sequence, and CTCF ChIP-seq to make cell-type-specific predictions, ChromaFold yields superior prediction performance when including CTCF ChIP-seq data as an input and comparable performance without. Finally, fine-tuning ChromaFold on paired scATAC-seq and Hi-C in a complex tissue enables deconvolution of chromatin interactions across cell subpopulations. ChromaFold thus achieves state-of-the-art prediction of 3D contact maps and regulatory interactions using scATAC-seq alone as input data, enabling accurate inference of cell-type-specific interactions in settings where 3C-based assays are infeasible.

Discovery of Competent Chromatin Regions in Human Embryonic Stem Cells.

The mechanisms underlying the ability of embryonic stem cells (ESCs) to rapidly activate lineage-specific genes during differentiation remain largely unknown. Through multiple CRISPR-activation screens, we discovered human ESCs have pre-established transcriptionally competent chromatin regions (CCRs) that support lineage-specific gene expression at levels comparable to differentiated cells. CCRs reside in the same topological domains as their target genes. They lack typical enhancer-associated histone modifications but show enriched occupancy of pluripotent transcription factors, DNA demethylation factors, and histone deacetylases. TET1 and QSER1 protect CCRs from excessive DNA methylation, while HDAC1 family members prevent premature activation. This "push and pull" feature resembles bivalent domains at developmental gene promoters but involves distinct molecular mechanisms. Our study provides new insights into pluripotency regulation and cellular plasticity in development and disease. One sentence summary: We report a class of distal regulatory regions distinct from enhancers that confer human embryonic stem cells with the competence to rapidly activate the expression of lineage-specific genes.

Dynamic network-guided CRISPRi screen identifies CTCF-loop-constrained nonlinear enhancer gene regulatory activity during cell state transitions.

Comprehensive enhancer discovery is challenging because most enhancers, especially those contributing to complex diseases, have weak effects on gene expression. Our gene regulatory network modeling identified that nonlinear enhancer gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Using human embryonic stem cell definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. We discovered a comprehensive set of enhancers for each of the core endoderm-specifying transcription factors. Many enhancers had strong effects mid-transition but weak effects post-transition, consistent with the nonlinear temporal responses to enhancer perturbation predicted by the modeling. Integrating three-dimensional genomic information, we were able to develop a CTCF-loop-constrained Interaction Activity model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and systematic enhancer discovery in both normal and pathological cell state transitions.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project.

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

Parallel genome-scale CRISPR screens distinguish pluripotency and self-renewal.

Pluripotent stem cells are defined by both the ability to unlimitedly self-renew and differentiate to any somatic cell lineage, but understanding the mechanisms that control stem cell fitness versus the pluripotent cell identity is challenging. We performed four parallel genome-scale CRISPR-Cas9 screens to investigate the interplay between these two aspects of pluripotency. Our comparative analyses led to the discovery of genes with distinct roles in pluripotency regulation, including many mitochondrial and metabolism regulators crucial for stem cell fitness, and chromatin regulators that control stem cell identity. We further discovered a core set of factors that control both stem cell fitness and pluripotency identity, including an interconnected network of chromatin factors that safeguard pluripotency. Our unbiased and systematic screening and comparative analyses disentangle two interconnected aspects of pluripotency, provide rich datasets for exploring pluripotent cell identity versus self-renewal, and offer a valuable model for categorizing gene function in broad biological contexts.

Dynamic network-guided CRISPRi screen reveals CTCF loop-constrained nonlinear enhancer-gene regulatory activity in cell state transitions.

Comprehensive enhancer discovery is challenging because most enhancers, especially those affected in complex diseases, have weak effects on gene expression. Our network modeling revealed that nonlinear enhancer-gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Utilizing hESC definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. The screen discovered a comprehensive set of enhancers (4 to 9 per locus) for each of the core endoderm lineage-specifying transcription factors, and many enhancers had strong effects mid-transition but weak effects post-transition. Through integrating enhancer activity measurements and three-dimensional enhancer-promoter interaction information, we were able to develop a CTCF loop-constrained Interaction Activity (CIA) model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and more comprehensive enhancer discovery in both normal and pathological cell state transitions.

MLL3 loss drives metastasis by promoting a hybrid epithelial-mesenchymal transition state.

Phenotypic plasticity associated with the hybrid epithelial-mesenchymal transition (EMT) is crucial to metastatic seeding and outgrowth. However, the mechanisms governing the hybrid EMT state remain poorly defined. Here we showed that deletion of the epigenetic regulator MLL3, a tumour suppressor frequently altered in human cancer, promoted the acquisition of hybrid EMT in breast cancer cells. Distinct from other EMT regulators that mediate only unidirectional changes, MLL3 loss enhanced responses to stimuli inducing EMT and mesenchymal-epithelial transition in epithelial and mesenchymal cells, respectively. Consequently, MLL3 loss greatly increased metastasis by enhancing metastatic colonization. Mechanistically, MLL3 loss led to increased IFNγ signalling, which contributed to the induction of hybrid EMT cells and enhanced metastatic capacity. Furthermore, BET inhibition effectively suppressed the growth of MLL3-mutant primary tumours and metastases. These results uncovered MLL3 mutation as a key driver of hybrid EMT and metastasis in breast cancer that could be targeted therapeutically.

Genetic Engineering of Immune Evasive Stem Cell-Derived Islets.

Genome editing has the potential to revolutionize many investigative and therapeutic strategies in biology and medicine. In the field of regenerative medicine, one of the leading applications of genome engineering technology is the generation of immune evasive pluripotent stem cell-derived somatic cells for transplantation. In particular, as more functional and therapeutically relevant human pluripotent stem cell-derived islets (SCDI) are produced in many labs and studied in clinical trials, there is keen interest in studying the immunogenicity of these cells and modulating allogeneic and autoimmune immune responses for therapeutic benefit. Significant experimental work has already suggested that elimination of Human Leukocytes Antigen (HLA) expression and overexpression of immunomodulatory genes can impact survival of a variety of pluripotent stem cell-derived somatic cell types. Limited work published to date focuses on stem cell-derived islets and work in a number of labs is ongoing. Rapid progress is occurring in the genome editing of human pluripotent stem cells and their progeny focused on evading destruction by the immune system in transplantation models, and while much research is still needed, there is no doubt the combined technologies of genome editing and stem cell therapy will profoundly impact transplantation medicine in the future.

CRISPR screening uncovers a central requirement for HHEX in pancreatic lineage commitment and plasticity restriction.

The pancreas and liver arise from a common pool of progenitors. However, the underlying mechanisms that drive their lineage diversification from the foregut endoderm are not fully understood. To tackle this question, we undertook a multifactorial approach that integrated human pluripotent-stem-cell-guided differentiation, genome-scale CRISPR-Cas9 screening, single-cell analysis, genomics and proteomics. We discovered that HHEX, a transcription factor (TF) widely recognized as a key regulator of liver development, acts as a gatekeeper of pancreatic lineage specification. HHEX deletion impaired pancreatic commitment and unleashed an unexpected degree of cellular plasticity towards the liver and duodenum fates. Mechanistically, HHEX cooperates with the pioneer TFs FOXA1, FOXA2 and GATA4, shared by both pancreas and liver differentiation programmes, to promote pancreas commitment, and this cooperation restrains the shared TFs from activating alternative lineages. These findings provide a generalizable model for how gatekeeper TFs like HHEX orchestrate lineage commitment and plasticity restriction in broad developmental contexts.

A hypomorphic mutation in Pold1 disrupts the coordination of embryo size expansion and morphogenesis during gastrulation.

Formation of a properly sized and patterned embryo during gastrulation requires a well-coordinated interplay between cell proliferation, lineage specification and tissue morphogenesis. Following transient physical or pharmacological manipulations of embryo size, pre-gastrulation mouse embryos show remarkable plasticity to recover and resume normal development. However, it remains unclear how mechanisms driving lineage specification and morphogenesis respond to defects in cell proliferation during and after gastrulation. Null mutations in DNA replication or cell-cycle-related genes frequently lead to cell-cycle arrest and reduced cell proliferation, resulting in developmental arrest before the onset of gastrulation; such early lethality precludes studies aiming to determine the impact of cell proliferation on lineage specification and morphogenesis during gastrulation. From an unbiased ENU mutagenesis screen, we discovered a mouse mutant, tiny siren (tyrn), that carries a hypomorphic mutation producing an aspartate to tyrosine (D939Y) substitution in Pold1, the catalytic subunit of DNA polymerase δ. Impaired cell proliferation in the tyrn mutant leaves anterior-posterior patterning unperturbed during gastrulation but results in reduced embryo size and severe morphogenetic defects. Our analyses show that the successful execution of morphogenetic events during gastrulation requires that lineage specification and the ordered production of differentiated cell types occur in concordance with embryonic growth.

Abnormal exocrine-endocrine cell cross-talk promotes ß-cell dysfunction and loss in MODY8.

MODY8 (maturity-onset diabetes of the young, type 8) is a dominantly inherited monogenic form of diabetes associated with mutations in the carboxyl ester lipase (CEL) gene expressed by pancreatic acinar cells. MODY8 patients develop childhood-onset exocrine pancreas dysfunction followed by diabetes during adulthood. However, it is unclear how CEL mutations cause diabetes. In the present study, we report the transfer of CEL proteins from acinar cells to ß-cells as a form of cross-talk between exocrine and endocrine cells. Human ß-cells show a relatively higher propensity for internalizing the mutant versus the wild-type CEL protein. After internalization, the mutant protein forms stable intracellular aggregates leading to ß-cell secretory dysfunction. Analysis of pancreas sections from a MODY8 patient reveals the presence of CEL protein in the few extant ß-cells. The present study provides compelling evidence for the mechanism by which a mutant gene expressed specifically in acinar cells promotes dysfunction and loss of ß-cells to cause diabetes.

Epigenome rewiring in human pluripotent stem cells.

The epigenome plays a crucial role in modulating the activity of regulatory elements, thereby orchestrating diverse transcriptional programs during embryonic development. Human (h)PSC stepwise differentiation provides an excellent platform for capturing dynamic epigenomic events during lineage transition in human development. Here we discuss how recent technological advances, from epigenomic mapping to targeted perturbation, are providing a more comprehensive appreciation of remodeling of the chromatin landscape during human development with implications for aberrant rewiring in disease. We predict that the continuous innovation of hPSC differentiation methods, epigenome mapping, and CRISPR (clustered regularly interspaced short palindromic repeats) perturbation technologies will allow researchers to build toward not only a comprehensive understanding of the epigenomic mechanisms governing development, but also a highly flexible way to model diseases with opportunities for translation.

Stage-specific regulation of DNA methylation by TET enzymes during human cardiac differentiation.

Changes in DNA methylation are associated with normal cardiogenesis, whereas altered methylation patterns can occur in congenital heart disease. Ten-eleven translocation (TET) enzymes oxidize 5-methylcytosine (5mC) and promote locus-specific DNA demethylation. Here, we characterize stage-specific methylation dynamics and the function of TETs during human cardiomyocyte differentiation. Human embryonic stem cells (hESCs) in which all three TET genes are inactivated fail to generate cardiomyocytes (CMs), with altered mesoderm patterning and defective cardiac progenitor specification. Genome-wide methylation analysis shows TET knockout causes promoter hypermethylation of genes encoding WNT inhibitors, leading to hyperactivated WNT signaling and defects in cardiac mesoderm patterning. TET activity is also needed to maintain hypomethylated status and expression of NKX2-5 for subsequent cardiac progenitor specification. Finally, loss of TETs causes a set of cardiac structural genes to fail to be demethylated at the cardiac progenitor stage. Our data demonstrate key roles for TET proteins in controlling methylation dynamics at sequential steps during human cardiac development.

Generation of human embryonic stem cell models to exploit the EWSR1-CREB fusion promiscuity as a common pathway of transformation in human tumors.

Chromosomal translocations constitute driver mutations in solid tumors and leukemias. The mechanisms of how related or even identical gene fusions drive the pathogenesis of various tumor types remain elusive. One remarkable example is the presence of EWSR1 fusions with CREB1 and ATF1, members of the CREB family of transcription factors, in a variety of sarcomas, carcinomas and mesotheliomas. To address this, we have developed in vitro models of oncogenic fusions, in particular, EWSR1-CREB1 and EWSR1-ATF1, in human embryonic stem (hES) cells, which are capable of multipotent differentiation, using CRISPR-Cas9 technology and HDR together with conditional fusion gene expression that allows investigation into the early steps of cellular transformation. We show that expression of EWSR1-CREB1/ATF1 fusion in hES cells recapitulates the core gene signatures, respectively, of angiomatoid fibrous histiocytoma (AFH) and gastrointestinal clear cell sarcoma (GI-CCS), although both fusions lead to cell lethality. Conversely, expression of the fusions in hES cells differentiated to mesenchymal progenitors is compatible with prolonged viability while maintaining the core gene signatures. Moreover, in the context of a mesenchymal lineage, the proliferation of cells expressing the EWSR1-CREB1 fusion is further extended by deletion of the tumor suppressor TP53. We expect the generation of isogenic lines carrying oncogenic fusions in various cell lineages to expand our general understanding of how those single genetic events drive tumorigenesis while providing valuable resources for drug discovery.

QSER1 protects DNA methylation valleys from de novo methylation.

DNA methylation is essential to mammalian development, and dysregulation can cause serious pathological conditions. Key enzymes responsible for deposition and removal of DNA methylation are known, but how they cooperate to regulate the methylation landscape remains a central question. Using a knockin DNA methylation reporter, we performed a genome-wide CRISPR-Cas9 screen in human embryonic stem cells to discover DNA methylation regulators. The top screen hit was an uncharacterized gene, QSER1, which proved to be a key guardian of bivalent promoters and poised enhancers of developmental genes, especially those residing in DNA methylation valleys (or canyons). We further demonstrate genetic and biochemical interactions of QSER1 and TET1, supporting their cooperation to safeguard transcriptional and developmental programs from DNMT3-mediated de novo methylation.

Enhancer Predictions and Genome-Wide Regulatory Circuits.

Spatiotemporal control of gene expression during development requires orchestrated activities of numerous enhancers, which are cis-regulatory DNA sequences that, when bound by transcription factors, support selective activation or repression of associated genes. Proper activation of enhancers is critical during embryonic development, adult tissue homeostasis, and regeneration, and inappropriate enhancer activity is often associated with pathological conditions such as cancer. Multiple consortia [e.g., the Encyclopedia of DNA Elements (ENCODE) Consortium and National Institutes of Health Roadmap Epigenomics Mapping Consortium] and independent investigators have mapped putative regulatory regions in a large number of cell types and tissues, but the sequence determinants of cell-specific enhancers are not yet fully understood. Machine learning approaches trained on large sets of these regulatory regions can identify core transcription factor binding sites and generate quantitative predictions of enhancer activity and the impact of sequence variants on activity. Here, we review these computational methods in the context of enhancer prediction and gene regulatory network models specifying cell fate.